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Abstract 

We analyze the fluctuations of the steady state profiles in the mod- 
ulated phase of the ABC model. For a system of L sites, the steady 
state profiles move on a microscopic time scale of order L 3 . The vari- 
ance of their displacement is computed in terms of the macroscopic 
steady state profiles by using fluctuating hydrodynamics and large de- 
viations. Our analytical prediction for this variance is confirmed by 
the results of numerical simulations. 

1 Introduction 

A surprising property of non-equilibrium systems in their steady state is that 
they may exhibit phase transitions in one dimension [H[2l[3l[ll[5l[6l[71[8l[9]. 
One of the simplest models for which such a phase transition occurs is the 
ABC model [El Q31 CEZ1 CGS El ESI QH 03 Q2] • In the ring geometry [TO1 Q31 
[T2l fl3l [T4] the model describes, on a ring of L sites, a lattice gas of three 
species of particles undergoing asymmetric exchanges between neighboring 
sites: each lattice site is occupied by a particle of type either A, or B or C 
and particles on neighboring sites exchange with the following rates 



AB ^± BA 

1 

BC ^± CB (1) 

1 

CA ^± AC 
1 
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where q < 1. The total numbers N A ,N B ,N C (with N A + N B + N c = 1) 
of particles of each species are conserved by the dynamics. Under these 
dynamics, the system reaches in the long time limit a steady state. When 
q = 1, all configurations are equally likely in the steady state. Therefore the 
steady state profiles are flat. As soon as q 7^ 1, the translational symmetry 
is broken and the steady state consists of three macroscopic domains corre- 
sponding to the three species \10\ lllj . The phase transition between these 
two regimes for the ring geometry has been studied in [121 113 1 [H I I19j. More 
recent works have also considered the ABC model on an open interval [15] 
and in its grand canonical version [TBI [T71 [18] . 
When the asymmetry q in ([JJ scales as 



exp 



(2) 



the dynamics is diffusive. This means that, on time scales of order L 2 , 
one can describe the system by density profiles {p A (x, r), p B (x, t), pc(x, r)} 
functions of a macroscopic coordinate < x < 1 and of a macroscopic time 
r which are related to the labels k of the lattice sites (1 < k < L) and to 
the microscopic time t by a diffusive scaling 

k = Lx ; t = L 2 r (3) 

In the ring geometry when one varies the asymmetry (3 in ([2]), one ob- 
serves a phase transition [10l CEH [12l [131 E] i n the steady state, from a flat 
phase where the density profiles {p A (x), pb(x), pc(x)} of the three species 
do not depend on the position x on the ring to a segregated phase where 
the profiles become modulated and therefore space dependent (on the ring 
all the density profiles are periodic functions of the macroscopic coordinate 
x, with period 1). As the dynamical rules in the ABC model are translation 
invariant, the continuous translational symmetry is broken in the segregated 
phase. Therefore, in the steady state, if one observes a set of three density 
profiles {p A (x), pb(x), pc(x)}, the same density profiles 

{p A {x-Z),p B {x-t),p c {x-i)} (4) 

translated by an arbitrary amount £ along the circle are also steady state 
profiles. For general densities of the three species, one does not know if the 
shape of the steady state profiles is unique. Throughout this paper, we will 
however assume that the uniqueness holds. 

In the thermodynamic limit, each set of the profiles @ labelled by a 
fixed value of £ are steady state profiles, in the modulated phase, with 
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an infinite life time and ergodicity is broken. On the other hand, for a 
large but finite system size L, the steady state is unique and the density 
profiles get translated by a time dependent random amount £ r due to the 
stochastic dynamics. Thus one observes in the steady state the following 
time dependence of the three profiles 

{PA(X - £t),Pb{x - ir),Pc{x ~ ir)} (5) 

where the phase £ T performs a stochastic motion along the ring. The goal 
of the present paper is to characterize the fluctuations of the phase £ r and 
to show that for a large system size its variance is given by 

Km I«a-(U 2 )-y (6) 

where L is the number of sites on the ring and the diffusion constant D is 
given in terms of the steady state density profiles (j4]) by 

D = 2 I dx (pApBPc + PBPCPa + PCPAPb) ^ 

[/ dx {pap'b- PaPb)] 2 

(One can easily check that this last expression is symmetric in p~a,Pb,PC 
using the fact that Pa + Pb + Pc = 1 which implies that Pap'b ~ P'aPb = 
PBp'c ~ P'bPc + Pb so ^at the squared integral in the denominator of ([7]) 
remains unchanged in any permutation of the three species.) 

The scaling of the variance ([6]) predicts a microscopic relaxation on a time 
t ~ L 3 which is exactly the same time scale found recently for this problem 
when the three species have equal densities [20J. Front fluctuations have 
been already analyzed (in full mathematical rigor) in |214 l22l [23] for one- 
dimensional non-conservative dynamics and in [23] for a system of coupled 
equations. The relation between front fluctuations and large deviations of 
forced interfaces has also been studied in [28J. 

The outline of the paper is as follows. In section 2, we briefly recall 
some known properties of the ABC model on a ring in the diffusive limit. 
In section 3, we present the calculation of the variance (J6] [7]) using fluc- 
tuating hydrodynamics. In section 4, we give an alternative derivation of 
the diffusion constant ([7]) based on large deviations in the spirit of [28]. In 
section 5, we show how our main result (JU [7J gets simplified in the equal 
density case (where the dynamics satisfies detailed balance) or when one 
approaches a second order phase transition. In section 6, we present the 
results of numerical simulations which agree well with the prediction ([6] [Jj . 
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2 The ABC model on a ring 



For the ABC model, on a ring of L sites, each site i is occupied by one 
particle of type A, B or C and the exchange rates between neighboring sites 
are given by ([1]). As the dynamics is ergodic and conserves the total number 
of particles of each species, the steady state properties depend only on the 
total densities r A ,r B ,rc of the three species which satisfy of course 

r A + r B + r c = 1 (8) 

and on the asymmetry q in the exchange rates ([I]). 

When this asymmetry q scales with the system size L as in (|2|), the dy- 
namics become diffusive and the density profiles pa(x,t), p B (x,r), pc(x,r) 
evolve in the infinite L limit according to [12] 



d T p A (x,r) = dlp A {x,T) + f3d x (p A {x,T)(p B (x,T) - pc(x,t))) 
d T p B (x,T) = dlp B (x,T) + pd x (p B {x,T)(pc{x,T) - pa(x,t))) (9) 
d T p c (x,T) = dlp c {x,T) + /3d x (p c (x,T)(p A (x,T) - p B (x,r))) 

on macroscospic scales ([3]). To describe the profiles on time scales much 
larger than L 2 , one would need to take into account stochastic corrections 
of order 1 /yL to these equations (see section [3|) and possibly higher order 
deterministic corrections in 1/L. 

Flat profiles {pa(x), p B (x), pc(x)} = {rA,r B ,rc} are always steady 
state solutions of these equations. It is however easy to see j!2j that they are 
linearly unstable when the asymmetry parameter j3 in (|2|) exceeds a certain 
critical value (3 C given by 

Pc = 7/2 10 

So when (3 > (3 C the steady state profiles are always modulated. 

In the equal density case (r^ = r B = rc = 1/3), it is known that (fTUI) 
gives the location f} c = 27r\/3 of the phase transition. For other densities, 
(jlOP is simply the value of /3 where flat profiles become linearly unstable, 
and so it coincides with the location of the phase transition between the flat 
phase and the modulated phase only when this phase transition is second 
order. One expects [12] in particular that a first order transition should 
occur at least when r A + r B +r c < 2(r A + r^ B + r c ) in which case (TL0j) is 
certainly not the location of the phase transition. 

The reason why the equal density case r A = r B = rc = 1/3 is best 
understood is that its dynamics satisfies detailed balance [10} [TT] and the 
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steady state measure is known. As a result the steady state profiles in 
the modulated phase do not move. For unequal densities {rA,rs,rc} 7^ 
{1/3,1/3,1/3}, the dynamics does not satisfy detailed balance and the 
steady state could move at a fixed velocity or even have a more complicated 
time behavior. We argue in Appendix I that the velocity should remain zero 
for some range of parameters in a neighborhood of the equal density case. 
Throughout this paper, we will assume that the steady state profiles do not 
move in which case they satisfy 

= d 2 x p A + (3d x (p A (p B -pc)) 

= d 2 x p B + pa x {p B {pc - pa)) (11) 

= dlpc + Pd x (pc(pa - Pb)) 

Expressions of these steady state profiles in terms of elliptic integrals 
have been given in [TU [131 E]- As each site is occupied by a particle A, B 
or C, one always has 

Pa(x,t) +pb(x,t) +pc(x,r) = 1 (12) 

Therefore one can restrict the analysis to the study of only two density 
profiles with Q replaced by 

d T p A = dip a + /3d x (p A (pA + 2p B - 1)) 

drPB = dips + /3d x (p B (l -2p A - pb)) (13) 

and (HU) by 

= dip A + pd x {pa{pa + 2p B - 1)) 

= dlp B + (3d x {p B (l-2pA-pB)) (14) 

When the system is large but finite, due to the stochastic nature of the 
microscopic dynamics, the density profiles pa(x,t), p B (x,r) (and pc(x,r) 
related to them by ()12p ) have no longer a deterministic evolution given by 
([9]). This deterministic evolution is replaced as in the macroscopic fluctu- 
ation theory |30|, [3"T| [32] by a probability distribution of observing macro- 
scopic density profiles {pa(x, t), pb{x, t)} during a macroscopic time inter- 
val t\ < t < T2 given by [33] 

Pro({pA,PB})~exp[-L min l[ TljT2] (p A , PbJa, 3b) ) (15) 
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where the large deviation functional is given by 

I [r u r 2 ](PA,PB,jAjB)= (16) 

T2 dr C dx aBB ^ A - qA ^ ~ 2ctab ( Jj4 ~ SvOO's _ 1b) + ^aa(Jb ~ Qb) 2 



with 



o"aa = o"aa(p) = 2pa(1 - Pa), 0"ab = <?ab{p) = -%PaPb 
°bb = ctbb(p) = 2p B (l ~ Pb) (17) 

and qA,QB given by 

QA = - f3p A (pA + 2/9 B - 1), gs = -^r- ~ PpbO- ~ %PA ~ Pb) (18) 

ax ax 

In (115j) the optimum is over the currents j A (x, t )j Jb(x, t) which satisfy the 
conservation laws 

dpA = _djA dfB_ = _djs_ , 1Q s 
dr dx ' dr dx 

An alternative way of representing the noise produced at the macroscopic 
scale by the stochastic dynamics of the microscopic model is to describe the 
evolution of the macroscopic profiles by fluctuating hydrodynamics [33]: the 
current profiles are then given by 

3A = QA + -j=tia(x,t), jb = QB + -j=t]b(x,t) (20) 

which together with (|19f) give the stochastic evolution of the macroscopic 
profiles. In (f20|) . i]A,rjB are Gaussian noises, white in time and in space, 
with the following correlations 

(va(x,t)) = (t] B (x,t)) = 

(rj A (x,T)r]A(x / ,T')) = cfaa{p(x,t)) 5(x - x) S(t - r') 
(r, A {x,T) VB (x',r')) = a AB {p(x,T)) 5{x - x') S(r - r') (21) 
{r] B (x,T)r] B (x',T')) = a BB (p(x,T)) 8(x - x') 5(-< 



(T-T 1 ) 



where the expressions of a AA , ct A b, obb in terms of the time dependent 
profiles p A (x,r), pb(x,t) are given in (fl"T)h 
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3 The fluctuating hydrodynamic approach 



For large L, in the steady state, the probability that the density profiles 
differ noticeably from the profiles p is exponentially small. The main idea 
to derive ([61 [7]) from the fluctuating hydrodynamic approach is to write for 
large L the density profiles as 

PA(x, t) ~ pA(x - £ T ) + -j=4a{x - £ T , T) 

Pb(x,t) ~ p B (x - f r ) + -^=cj) B {x - £ T ,r) (22) 

where pA and p B are the steady state profiles of the infinite system, 4>a 
and (f>B are the fluctuating parts of these profiles due to the noise and £ T is 
the cumulative translation of the profiles due to the noise. (pA an d do 
not grow with time. The approach we then follow is inspired by the works 
[2~H |24"1 |2~5| [26l [27] on the effect of noise on the position of fronts. 

In order to manipulate conveniently these quantities, it is useful to in- 
troduce the vectors 

-/ \ ( Pa{x) \ if \ ( <Pa(x,t) \ f \ ( Va(x,t) 

p{x) =[_}'] ; <p(x, r) = ; ' ; rj{x, r = ) ' 



and the covariance matrix of the noise (|171 \21\i 



(23) 



= ( ^rill ) ( 24 ) 

V o-ab(p(^)) <tbb\P\x)) J 

Let us also define the operator £, acting on such space dependent vectors, 
which is obtained by linearizing the macroscopic equations (1131) around the 
steady state profiles p 



C 



ipA{x) \ ( 9 x^a(x) +Pd x U2pA(x) + 2p B (x) - 1)^a{x) + 2pA(x)ip B (x) 
M*) )~y 8 2 x i{] B (x) + /3d x ((l - 2p A (x) - 2p B (x))^ B (x) - 2p B (x)7p A (x) 

(25) 

It is easy to check using (fl4l) that 

C p'(x) = (26) 

as expected from the translation invariance @. 

As the dynamics conserve the total number of particles of each species, 
the fluctuating parts 4>a{x,t) and 4> b (x,t) satisfy 

4>a(x,t)cIx = / 4> B (x, r)dx = (27) 
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Let us assume that p'{x) is the only eigenvector of C with a zero eigenvalue 
in the space of functions which satisfy (|27p and that all the other right 
eigenvectors tp^ of C (for a > 1) have eigenvalues X a with negative real 
part 

C ^ a \x) = X a ^ a \x) (28) 

In the following we will also use the scalar product between a left vector 

tpA(x) 



x( x ) — (xa(x), Xb(x)) and a right vector i/j(x) 



4) B {x) 



<x|V>> 



i r 



XiW^W + Xb(x)iI) B (x) 



dx 



It is easy to check, using ([25]) . integrations by parts and JT 
right vector i/j(x) 

(x {0) \c^) = o 

where the vector x^°\ x ) is given by 

x (°)( x ) = (p B (x),-pA(x)) 



(29) 

that for any 
(30) 

(31) 



So (|3ip is a left eigenvector of C with zero eigenvalue. It follows that for all 
the eigenvectors ijj( a \x) with a non-zero eigenvalue X a one has 



< 



(32) 



If we put the expressions (|22p into the equations of evolution of the 
density profiles (fTBl [T9|) one gets at leading order in 1/vT 

1 11 

- (drCr)p(x) + — d T (p(x, r) = — Ccj)- —d x 7](x, t) (33) 

From (|33p . one can see already that the variance of £ T scales as 1/L as 
claimed in ([6]). This L dependence is consistent with the L 3 dependence of 
the relaxation time found in [20] in the case of equal densities (r^ = r B = 
rc = 1/3) (see section 6 and figure 1 below). 

Under the assumption that one can decompose d x r](x, r) on the eigen- 
vectors of C and that p'(x) is the only eigenvector with zero eigenvalue 
satisfying (I27p . one can write 



d x r,(x,r) = c (T)p'(x) + ]T c a {r)^ a \. 
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This together with (|33|) leads to 



0( x ,r) = -^V (Q) (x) / ^(rOe^^'^r' 

Q>1 

and to 

d^r = ^ (34) 

Doing so, we have imposed that (j) has a zero component on p'(x) (the com- 
ponent on p'(x) of the r.h.s. of (|33p can always be absorbed as a translation 
which contributes to £ T ). 

As <^(x) has only components on the right eigenvectors of C with non 
zero eigenvalue, one can determine cq(t) by taking the scalar product of (|33|) 
with the left eigenvector ( x ) an d by using the orthogonality property (|30| 
[32]) one gets 

o £ J_ : x (0) (^)l^r?(x,r) > 
which after an integration by parts becomes 

At ~ VI <yW(x)\p{x)> [6b) 



This together with (|17|21l24p leads to the following expression for the diffu- 
sion constant D defined in 

= <d xX {0 \x)\^(x)\d xX ^(x)> 

<X i0) {x)\p'(x)> 2 [ ' 

where we have used the fact (|21I24|) that S(x) is the covariance matrix of 
the noise tj(x,t). From ([24"1 [3Tj) one deduces 

D = I dx(2p A (l - Pa)Pb + ^PaPbP'aP'b + - Pb)p'a) 

[I dx (pap' b - p'aPb)] 2 

which can be rewritten as ((7]) using (fT2j) . 

Although (|22|) should remain valid for times t > I, we have neglected in 
the derivative (|33l36p of £ r corrections of order 1/L. Integrating (|36|) gives 
£ r ~ y^x wmcn might compete with a possible drift correction of order . 
Therefore the value of £ T obtained by integrating (|36p should only be valid 
for r < L. 
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4 The large deviation approach 



In the large deviation approach we compute the probability (|151 [16]) of ob- 
serving profiles close to the steady state profiles Pa(x), Pb(x) moving for a 
long time at a given velocity v around the circle. 

For small v one expects these profiles to be of the form 

Pa(x, t) = pa(x - vt) + v iPa(x - vt) + 0(v 2 ) 

Pb(x, t) = Pb{x — vt) +v i/jb(x — vt) + 0(v 2 ) (38) 

where tpA and tps are the deformations of the profiles due to the velocity v. 
we have to look for the deformations ipA, ^B which maximize the probability 
of such moving profiles. To produce such time dependent density profiles, 
the currents (I19p should be of the form (at first order in v) 

j A (x, t) = I A + v i A + v p A (x - vt) + 0(v 2 ) 

j B (x, t) = I B + v i B + v p B (x - vt) + 0(v 2 ) (39) 

where Ia,Ib are the steady state currents, which satisfy 
I A = -d x pA - Ppa(pa + 2p~B - 1) 

Ib = -d x p B - PPb{1 - 2p A - Pb) (40) 

lAi^B are additional currents and the terms vpA(x—vt),vpB(x—vt) are there 
to insure that the conservation laws (|19p are satisfied. For small velocity 
v, the cost of the large deviation functional (I16p evaluated for the moving 
profiles (|38p can be written as 

c y 2 f 1 dx Pb(1 - PB)(iA + PA- <Pa) 2 Pa{1 ~ PA)(iB + Pb ~ 0b) 2 
J 4paPb(1 - Pa- Pb) 4paPb{1 - Pa- Pb) 

| 2 PAPB(iA + PA- 4>A)(iB +PB ~ 4>b) ^ 
4pAPfi(l - Pa- Pb) 

where 

4>(x) = Mtp(x) (42) 
-d x ^ A {x) - p((2p A (x) + 2p B (x) - 1)4>a{x) + 2p A (x)i> B {x) 
-d x i/j B (x) - P[ (1 - 2p A (x) - 2p B (x))ip B (x) - 2p B {x)ij} A (x) 



is obtained by linearizing (|18j) . One can notice from (|25p that for any right 
vector ip(x) 

d x (Mtfj(x)) = -£4>{x) (43) 
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This implies that (see (J31I32D ) that d x x {0) (x) = (p' B (x),-p' A (x)) is a left 
eigenvector of Ai with a zero eigenvalue and therefore that for any right 
vector ip(x), one has 

< d xX m {x)\M\ip(x) >= (44) 
Using the scalar product (J29]) one can rewrite (|4ip as 



v 2 

C = — <i + p(x) - <f>(x) \Z-\p(x))\ i + p(x) - <t>{x) > (45) 

where the inverse of the matrix Y>(p(x)) defined in (|24p is 

1 / Pb(x)(1 - Pb{x)) p A (x)p B (x) 



S~ (p(x)) = 

2p A (x)p B {x)p c {x) V pa(x)pb{x) p A (x)(l - p A (x)) 

(46) 

with pc{x) = 1 — Pa{x) — pb{x). We have now to optimize this expression 
over the currents i and the deformations ip which appear in ()45|) through 
(j4~2j) . To do so let us consider the right vector d x (T l (p(x))d x x^°\x)) . As it is 
a derivative, this vector has a zero average (j27|) on the circle and we assume 
that it can be decomposed on the basis of the eigenvectors of C 

d x (^{p(x))d xX {0) (x)) = b p'(x) + Yl h M*) = bop(x) + C${x) 

where ip(x) = ip a (x) or equivalently 

p\x) = ad x (Z(p(x))d xX {0 \x)) + C^{x) (47) 

(with a = 1/bo and i/j = —tp/bo). Integrating over x using (H3]l one gets 

p(x) = aZ(p(x))d x x {0) (x) - Mip(x) + i (48) 

where i is just the constant of integration. Then replacing p(x) by (|48p into 
U@5\i , one gets using the orthogonality property (|4"4"|) 

v 2 r 

C = — y< i + i - M(tp + 1 {p(x))\i + i - M(ip + $)> 

+a 2 < d x x i0) (x)\Z(p(x))\d xX {0) (x) >] (49) 

As the matrix S _1 (see dM|)) is positive definite, it is easy to see that the 
optimal choice for i and ip is i = — i and ip = — ip. By multiplying (|48p by 
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the left vector d x x^°\x) and by using the orthogonality property (|4"4"]). one 
gets for the projection a in (|48p 

<^x (0) (^)Ip(^) > 



< d x x^Kx)\E(p{x))\d xX ^Hx) > 

Therefore replacing i,tp,a by their expressions into (|49p gives for the prob- 
ability of seeing moving steady state profiles at a small velocity v during 
some macroscopic time r 

L v\ <d xX (0) (x)\p(x)>' 



Pro(-u) ~ exp ( — LminC ) = exp 

ip,i 



2 <d xX ^(x)\j:(p(x))\d x x^(x)> 



In appendix II, we give another derivation of (|50p in the reversible case 
ta = i~b = rc = 1/3. 

Remark 4.1 Instead of the scalar product (|29p . it could be also convenient 
(see for example \29j) to consider the modified scalar product 

(/U)-i,s = (V- 1 /|S- 1 (p)|V- 1 5 > (51) 

V^f(x) = [ X duf(u) + K, (52) 
J o 



where V 1 f is defined for vectors f(x) = with zero mean (f27|) by 



and the constant K is fixed such that (S(p) 1 V 1 /) = 0. 
The large deviation cost (|45p can be rewritten 



v 2 



C = —{d x p - d x (f> | d x p - d x (f))_ ix 

The optimal value of i in (|45|) is implicitly chosen by the condition (|52p . 
From the identity (|43p . one finally gets 

c = — (p - Cip | p - A/y_ 1)S 

27ie latter expression provides another way to understand the decomposition 
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The large deviation result (|50p was obtained in the limit where r and 
L go to oo first and then v tends to 0. If we rewrite (|50p in terms of the 
displacement £ r = vt, we get 



Pro(£ r ) ~ exp 



L e <d xX (o) (x)\p(x)> 2 

2t < d xX ^(x)\^(p(x))\d xX ^Kx) > 



(53) 



Formally, this leads to the same expression for the diffusion constant 
([6]) as (|3T|) . The difficult to recover rigorously the diffusion coefficient is 
that (|53p has been obtained for £ r of order r while to extract the diffusion 
constant from ([53]) one would need £ T to be of order y^j, that is for v of 

order \J~^j i - The fact that (|53p does give the right diffusion constant means 
that some kind of exchange of the order of limits r, L — > oo and v — > 
should hold. 

As we will see in the simulations of section [6] (Figure 2) , £ r may have a 
drift of order t/L which is not predicted by (|53l) . 

5 Two particular cases 

We now briefly discuss two cases where the expression ([7]) simplifies. 
5.1 The equal density case 

When the densities of the 3 species are equal (r^ = rs = tq = 1/3), the 
dynamics satisfy detailed balance. Therefore the currents I a and Ib in (|40p 
vanish and the steady state profiles satisfy 

d x PA = P{Pa{Pc ~ Pb)) 

d x PB = /3{pb(pa - Pc)) (54) 
d x pc = (3{pc{pb - pa)) 

This implies that d x {pApBPc) = and 

Pa(x) Pb(x) pc{x) = r (55) 

where T is a constant as already noted in |15] , 

By replacing p' c (x) 2 by /3p' c (x)pc(x)(pB{x) — Pa(x)) and doing similar 
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changes for p' A (x) 2 and p' B {x) 2 in the numerator of ([7|) one gets 

2 J dx (pApBPc + PbPcPa + PcPap'b) 

= 2/?r J dx[p' c (x)(p B (x) - pa(x)) + p' A (x)(p c (x) - p B {x)) + p' B (x)(p A (x) - pc(x))} 
= 6/?r J dx[p A (x)p B (x) - p B (x)p' A (x)} 
Therefore the ratio ([7]) becomes 
D 



6pr 



f dx[p A (x)p' B {x) - p B (x)p' A (x)} 
one can further transform the denominator using (|54p 

J dx[p A p B - PbPa] = P j dx{p A pB(l ~ 3pc)) = f3T 

Finally the expression ([7|) becomes 

dx 



dx 
Pc{x) 



D = 6 



Pc{x) 



(56) 



(in the equal density case (|56|) would of course remain unchanged by replac- 
ing pc by pA or p B ). 

5.2 Close to the second order phase transition 

When the phase transition of the ABC model on the ring is second order, and 
(3 approaches /3 C from above, the steady state profiles become sine functions 



p A (x) ~ r A + * 
Pb(x) c± r B + ^ 

pc(x) ~ r c + * 



r A eKT(*-*o) + c . c . 



- ^ - r B - ~ 2(r A + r\ + r, 2 



A'B'Ci 



tb ~ r A ~ r c + i\/ 1 - 2(r^ + r| + r 2 . 



2^ 



e 2«r(x-x ) + cc 



14 



where 



2(r 2 A + r% + r 2 c ) - 4(r 3 A + r% + r%) 



Then the integrals which appear in ([7]) become for (3 close to f3 c 
2 J dx (papbPc + PbPcP'a + PcPap'b) — ^ 2 r A r B r c * 



and 



So that the diffusion constant D diverges as (ft — (3 C )~ 1 
SrATBrc 



D 



(1 - 2r\ - 2r% - 2r 2 c )t> 2 

_ 6r A r B rc(r 2 A + r| + r 2 c - 2(r\ + r% + r^)) /3 C 



(l_ 2 ( r 2 +r | + r 2))3 

6 Simulations 

From (|5|) one can deduce the long time decay of the autocorrelation function: 
if Sk(t) = A, B or C is the type of particle on site k at time t, then in the 
steady state, for times \t — 1'\ = L 2 t » L 2 

Prob(S fe (i) = S k (t')) = (58) 



dxp A (x)pA(x - £t) + pb(x)pb(x - £ r ) + pc(x)p c {x - i T ) 



where the expectation Eg T is over the Gaussian variable £ T . Assuming that 
£ T has no drift and behaves as a Brownian motion with diffusion coefficient 
([6j [7J) then one gets in microscopic units (|3j) 



Prob(5 fe (t) = S k (t')) - (r\ + r% + r 2 ,) 

= ^(|a n | 2 + |6 n | 2 + |c n | 2 )e- 2 " 2 ^/ i 
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In (|59|) the a n are the Fourier coefficients of the steady state density profiles 



Pa{x) = r A + 



(60) 



(with similar definitions for b n and c n ) and r^, r#, re are the global densities 
of the three species. Should £ T have a drift of order 1/L (as measured in 
the non equilibrium case), then additional oscillations would occur in the 
autocorrelation function ([59]) • 

In the equal density case (r^ = re = rc = 1/3), we have performed 
Monte Carlo simulations of the ABC model to measure the autocorrelation 
function predicted by (|59p . At each time step we also measured the first 
Fourier component of the species A 



and the autocorrelation {cos(V(t') — V(t))) of the phase V(t). The phase 
V(t) is related to the displacement ^ T by V(tL 2 ) = 2ir ^ T where the equality 
is modulo 2n. Our results are shown in a semilogarithmic scale in Figure 
1 for three different system sizes L = 60, 120, 240 when j3 = 15. The data 
for the three different sizes overlap, indicating that the right time scale is 
L 3 . There is also an excellent agreement with the prediction (|7|), which in 
this case would give D ~ .94 estimated by an independent calculation of the 
profiles obtained by looking at the solution of the deterministic equations 



In attempts to check that the steady state profiles have no velocity, as 
claimed in the appendix I, we performed simulations of the ABC model for 
3 system sizes (L = 60,120,240) at densities ta = 1/3, rs = 1/2, rc = 1/6 
and /3 = 17. The position £ T of the profiles was measured by following 
the phase V(t) = 2tt^ t of the first Fourier mode of the density profile of the 
species A f)61 1) . Our results in Figure 2 (left) show that for these densities the 
displacement £ r has a linear time dependence and therefore the profiles have 
a non-zero velocity, but this velocity decreases with the system size. When 
we multiply the displacement by the system size L as in Figure 2 (right), we 
see that the three sets of data overlap indicating that the velocity scales as 
1/L. We observed this 1/L decay of the velocity for other choices of the den- 
sities, but the decay looked slower when we approached the phase transition 
line. We did not succeed in computing this velocity vl{ta-, rs, rc) but we 
think that it should be an antisymmetric function of the three densities and 



L 




(61) 



Id- 



le 



Prob (S k (t)=S k (f)) 




It'-tl/L 3 



Figure 1: The exponential decay of the autocorrelation Prob(iSjfe(t) = Sk(t')) (lower 
curves) and of (cos(V(t') — V(t))) versus \t — t'\/L 3 in the equal density case, for 
j3 = 15, and for 3 system sizes: L = 60, 120, 240. The thin line indicates the slope 
predicted by our main result ([7]). 

therefore should vanish whenever two densities are equal. Understanding 
this velocity of order 1/L would require to rewrite the fluctuating hydrody- 
namics equations (|33p to a higher order in 1/L or to calculate corrections 
of order 1 in the large deviations (|15p . For Ginzburg-Landau equation with 
white noise, it has been shown in [23] that the interplay between the non- 
linearity and the noise may induce such a small drift. It would be interesting 
to see if an analogous approach could be used for the ABC model. 

7 Conclusion 

In this paper, we have studied the fluctuations of the steady state profiles of 
the ABC model in the segregated phase. We have shown that the position of 
the density profiles fluctuates on a microscopic time scale of order L 3 and its 
variance has been computed in terms of the steady state profiles ([7]). The 
result has been obtained by using fluctuating hydrodynamics and a large 
deviation approach. It has been confirmed by numerical simulations (see 
section [6]) . 

Our result relies on the assumptions that the solutions of the macro- 
scopic steady state equations (fl~3j) are unique (up to translations) and do 
not move. This claim is correct at equal densities but it would be important 
to determine the range of parameters for which it remains valid. Both meth- 
ods we used (linearized fluctuating hydrodynamics and large deviations) are 
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Figure 2: On the left, the linear displacement of £ T for three system sizes L = 
60, 120,240 at densities r A = 1/3, r B = 1/2, r G = 1/6 and /3 = 17. On the right, 
the same data are depicted after rescaling by 1/L. 

valid for macroscopic time scales and as the fluctuations of the steady state 
occur on a longer time scale we had to extrapolate the results given by these 
methods. Several questions remain in order to understand the validity of 
this extrapolation as well as the intriguing small drift of order 1/L observed 
in the numerical simulations (see figure 2). 

It is our pleasure to dedicate this work to Professor Cyril Domb on the 
occasion of his 90th birthday. 
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Appendix I 

In this appendix, we argue that for densities rA,rs,rc close to the equal 
density case, the steady state profile solution of (fT3l) do not move. 

Let us assume that the steady state profiles of the macroscopic equations 
(fl~3j) have the form p(x — vt) = {p~a(x — vt),pb(x — vt)}. When = rs = 
1/3, the dynamics are reversible. The invariant measure is known and the 
steady state profiles have no velocity and are solutions of (fT4"|) . We are 
going to check that this property holds for general mean densities provided 
smoothness assumptions are satisfied. 
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Suppose that for the mean densities {rA, r B}, the steady state profile 
p = (^j has a velocity v and is solution of f)13|) . For a small shift of the 

densities {ta + Sa, tb+^b}, we assume that the steady state profile p = (^) 
is obtained as a perturbation of p 

Pa(x,t) = p A (x - v 5 t) + ijj A (x - v s t) ^ 
p B (x, t) = p B {x -v s t) + ip B (x - v s t) 

where v s — v and ip = Q^) are of order 5. As p is also a solution of (|13p . 
one gets at the first order in 5 

— (vs — v)p — vtp = Cil> (63) 

where the linearized operator C has been introduced in (|25p . By construction 
p' is a right eigenvector of C = C + vd x with zero eigenvalue. Making the 
assumption that there is a corresponding left eigenvector \^ of C with zero 
eigenvalue. Then 

(v s -v) <X {0) \P >=0 

Thus the derivative of the velocity with respect to 8 is zero since < x^\p' 
0. If one can interpolate by a smooth path {r^(s), rB(s)}o<s<i from {1/3, 1/3} 
to {va^b} and iterate the same argument along this path, then the veloc- 
ity of the steady state remains equal to zero. This argument relies on the 
assumption that the profiles and the linearized operator behave smoothly 
with respect to small density shifts. 

Let us finally add that we solved numerically the equations (JX3|) for 
several choices of the global densities ta^b and in all cases the solutions 
converged to non moving steady states. 



Appendix II 

For reversible dynamics (ta = ^b = 1/3), the dynamical large deviations 
(|15p can be related to the steady state large deviations which are known 
explicitly [12]. In the steady state, the probability of observing macroscopic 
density profiles pa(x), Pb{x) has the following large L dependence 

Pto(pa(x),pb(x)) ~ exp ( - LF[p A {x), Pb{x)} ) (64) 



19 



with 

F[p A {x),p B {x)]= (65) 

k+ / dx [p A (x) In p A (x) + p B (x) In p B (x) + pc (x) In pc (x)] 
Jo 

+/3 dx dz z[p B (x)p c (x + z) + pc(x)pa(x + z) + pa(x)pb(x + z)] 
Jo Jo 

where pc( x ) = 1 — Pa{x) — Pb(x) and k is a normalisation constant. As 
noted in the appendix of |33| . one can rewrite the evolution (fT3j) of the 
density p = (^) in terms of the large deviation functional T 

d T p(x,T) = d x (v(p(x,T))d x F(x,T)") (66) 



where F = ( F A ) is the vector defined by 



F A (x, r) = -^- (p(x, t)) , F B (x, t) = -^- (p(x, r)) 
bp a op B 

and X(/o) was introduced in (|24p . In particular the typical currents q = (^) 
defined in (fT8|) satisfy 

g (x,T) = -^(p(x,r))F(x,r) (67) 

We remark that for p = p then F = and q = 0. 

The previous identities allow us to rewrite the functional (|16p for the 
joint deviation of the density p = (^ A ) and the current j = (j A ) as 



1 j-T . 

Z[o,Ti(p,j) = 2j Q dr V ~ q S ^ 



Expanding the functional one obtains 

The cross product simplifies by using (I67D and dtp = —d x j 

2 / T dr(g | S" 1 ^) | j> = - f dr{H{p)d x F | JT^p) | j) 
j o Jo 

= -J dr{d x F\j) = -J dr(^\d t p\ = -T(p T )+T(p ) 
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Thus one has 

%T](p,i) = \{HPT) - HPo)) + \ £dr (j | E" 1 ^) | j) + (q \ ^(p) \ q) 



The expression (|68p of the large deviation functional shows that for small 
perturbations around the steady state profiles, the contributions of the cur- 
rent deviations and the profile deformations decouple. Observing a travelling 
wave p(x — vt) of the form (f38j) imposes a current j(x, r) = vp(x — vt) + vi 
at the first order in v (see ([39]) ). When T diverges, (j68|) becomes at the 
second order in v 

1 v 2 _ 

C = ^fi^ -f X [o,T]{p,j) = y (p + i | ^(p)^ 1 | p + i) + (g | S(p)^ 1 | g) 

where g was defined in (I67p in terms of As is a non negative quadratic 
form, we see that the lowest cost is achieved when the profile p is not mod- 
ified, i.e. when tp = 0. Thus it remains to optimize over the constant i and 
the large deviation cost at the second order in v is given by 

C= lim ±Z [0>T] (pJ) = ^(p + i\Z(p)- 1 \p + i) (69) 

with 

i = -(S(p)- 1 >- 1 <S(p)- 1 p) (70) 

One can easily check that (|69|) coincides with the expression obtained in 
([50]) when ta = tb = tq = 1/3. 
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